Quantum discord in spin systems with dipole-dipole interaction 
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The behavior of total purely quantum correlation (discord) in dimers consisting of dipolar-coupled 
spins 1/2 is studied. We found that the discord Q = at the absolute zero temperature. With 
increasing the temperature T, at first the quantum correlations in the system increase, smoothly 
reach the maximum, and then turn again into zero according to the asymptotic law . It is also 
shown that in absence of external magnetic field B, the classical correlation C at T — !> is, vice versa, 
maximal. Our calculations predict that in crystalline gypsum CaS04 ■ 2H2O the value of natural 
{B — 0) quantum discord between nuclear spins of hydrogen atoms is maximal at the temperature 
of 0.644 ^K, and for 1,2-dichloroethane H2CIC — CH2CI the discord achieves the largest value at 
T = 0.517 fiK. In both cases, the discord equals Q « 0.083 bit/dimer what is 8.3% of its upper limit 
in two-qubit systems. We estimate also that for gypsum at room temperature Q ~ 10~^* bit/dimer, 
and for 1,2-dichloroethane at T = 90 K the discord is Q ~ 10~^^ bit per a dimer. 

PACS numbers: 03.65. Ud, 03.67.Mn, 75.10.Jm, 75.50.Xx 



> 

m 
o 



X 



I. INTRODUCTION 

The notion quantum discord has been intro- 
duced in the quantum information theory by 
Zurek^. Discord was regarded as "a measure 
of a violation of classicality of a joint state of 
two quantum subsystems"—. However later its 
initial definition undergone some changes. 

In 2001, Henderson and Vedral and then in- 
dependent of them OUivier and Zurek in their 
papers^i^ (see also^i^) performed analysis of ev- 
ery possible correlations J in a bipartite system 
and suggested the ways to extract from them, 
on the one hand, the purely classical part C 
and, on the other hand, only the quantum con- 
tribution Q. The quantum excess of correla- 
tions, Q ^ I — C, has been called "discord" — 
in the modern understanding of this tern*^. 

The authors of above papers established also 
that the quantum correlation can be non-zero 
even in separable (but mixed) states. In other 
words, quantum correlations are not exhausted 
by entanglement {E). Entanglement, which 
can relate the different parts of a system even 
when there are no interactions between these 
parts (the Einstein-Podolsky- Rosen effect), is 
only a special kind of quantum correlations. 
Since the 80-90s of past century the entangle- 
ment was considered as a fundamental resource 
for quantum information processing, teleporta- 
tion, cryptography, metrology, and other tasks 
in quantum technology^^^'^. 

It is remarkable that quantum discord can 
also lead to a speedup over classical compu- 
tation and lead even without containing much 
entanglementiiii^. This important property 
of discord evoked extremely great interest to 



the new kind of correlations. One discovered 
also that discord can detect the quantum phase 
transitions^-ii^. Moreover, it has been shown 
that in contrast to entanglement and thermo- 
dynamical quantities, the discord makes it pos- 
sible to catch the approach of quantum phase 
transitions even at finite temperatures^^. Other 
important features of discord have been also 
noted. By this, a surprising fact turned out: 
"Almost all quantum states have nonclassical 
correlations"-^. Achieved up to now results on 
the theory and applications of quantum discord 
are given in the recent leviewsi^^. 

The goal of this paper is to study the be- 
havior of discord in spin systems with dipo- 
lar couplings. There is a large number of sub- 
stances magnetic interactions in which have the 
dipole-dipole character, and exchange and indi- 
rect ones are weak enough. By this, the spins 
both electrons and nuclears can serve as ele- 
mentary magnetic moments. The class of dipo- 
lar magnets with electron spins includes, for ex- 
ample, the Tutton salts, alumsi^, and numer- 
ous salts of rare earth elements^S. The typi- 
cal temperatures at which the effects of their 
dipole interactions show themselves lie in the 
millikelvin region. However the spin-lattice in- 
teractions in electron paramagnets are strong 
that leads to short relaxation times too. Nu- 
clear spins, of course, have an indubitable asset: 
the spin-lattice relaxation times for them can 
achieve minutes and hours. Using the available 
NMR data for the two classical examples, gyp- 
sum CaS04 • 2H202i and 1,2-dichloroethane 
HzClC-CHzClS^ (see also, e.g.^^) which 
contain the sufficiently isolated pairs of dipolar- 
coupled nuclear hydrogen spins, we estimate 
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the discord between these spins. 

In following sections of the paper we give 
definitions for the different correlations, formu- 
late the model, calculate the classical and quan- 
tum information correlations in it, estimate the 
discord for materials with spin-nuclear dimers, 
and, lastly, briefly summarize the results ob- 
tained. 



II. CLASSICAL AND QUANTUM 
CORRELATIONS 



In statistical theory, a degree of relationship 
(correlation) between two random variables x 
and y with the joint probability distribution 
function p{x, y) is often measured by covaria- 
tions or by Pearson's correlation coefhcient 



R = 



{x - x){y - y) 



(1) 



Here the bar denotes the average over the p rob- 
ability distribution and Di = (x — xY and 
D2 = [y — yY are dispersions respectively for x 
and y. Notice that in the mathematical statis- 
tics other types of correlation coefficients are 
used also. One should emphasize that the con- 
dition i? = does not imply, generally speak- 
ing, the independence of random variables, i. e., 
that p{x,y) = pi{x)p2{y)r^- 

In the classical information theory 
(e. g.,— 1^), one uses the notion of mutual 
information 

I{X : Y) = H{X) + H{Y) - H{X, Y) > (2) 

between two objects X Y. Here H{X), H{Y), 
and H{X, Y) are the Shannon entropies 

i/(X)--^Pi(x)logpi(x), 

X 

H{Y) = -Y,P2{y)\ogp2{y). 

V 

H(X,Y)^-Y,P{^.y)^ogp{x,y), (3) 

'x,y 

where pi{x) = Y^yPi^^y)^ P^iy) = J^xPi^^v)- 
(A choose of logarithm base defines the infor- 
mation measure unities: bits, nats, dits, hart- 
ley.) It is remarkable that now the equality 
/ = is a necessary and sufficient condition 
for the independence of X and Y. This prop- 
erty allows to use the mutual information as a 
measure of information correlation between the 
systems X and . 



Taking into account the Bayes rule one can 
rewrite the right part of Eq. ([2]) in Shannon's 
nonsymmetric form 



where 



I'iX ■.Y)^H{X)-HiX\Y), (4) 



H{X\Y)^H{X,Y)~H{Y) (5) 



is the condition entropy. In the classical case, 

/' EEE /. 

In the quantum information theorySiS2^ the 
equation ([2|) is replaced by the new definition 

I{A : B) = S{pa) + S{pb) - S{pab), (6) 

which serves as a measure of quantum mutual 
information between the two subsystems A and 
B composing together the joint system AB = 
A\J B. In Eq. (jS]), pab is the density matrix of 
the joint system AB^ pA and pB are the reduced 
density matrices, respectively, for subsystems 
A and B, and S{p) {p = {pA, Pb, Pab}) repre- 
sents the von Neumann entropy 



Sip) = -Tvp\ogp. 



(7) 



It is important that / = is the necessary and 
sufficient condition for the factorization pab = 
Pa ® Pb, what means, of course, the absolute 
independence (non-correlativety) of A and B in 
the product state. Therefore, in the quantum 
information theory, one takes the mutual infor- 
mation to measure the total (both classical and 
quantum) correlations between the two subsys- 
tems of bipartite quantum system. 

On the other hand, measurements performed 
on one system, in general, influence on the 
quantum state of another system (see, for 
instance, ^°). Postulating that the total classi- 
cal part of correlations is the maximal amount 
of information about one subsystem, say A, 
that can be extracted by performing a measure- 
ment on the other subsystem B, Henderson and 
VedraP suggested to take as a measure of clas- 
sical correlation the quantity 



C{pab) = max{S{pA) - Ylp^S{p'A)}. 



(8) 



Here {Bi} is a complete set of measurements 
on the subsystem B, 

p\ = Ttb{B.,pabB+)/Ttab{B^pabB+) (9) 

is the remaining state of A after obtaining the 
outcome i on B, and 



Pi = TvABiBiPABB:^ 



(10) 
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is the probability to detect the result i. 

OUivier and Zurek^,, on the contrary, focused 
their attention on an extraction of quantum 
correlations. Further analysis of the measure- 
ments led to the generalization of expression (|4]) 
for the quantum case, 



I'{A:B)^S{pa)~Y.P^^(p'^ 



(11) 



(It is obvious that the right hand of this equal- 
ity is a non-optimized classical correlation of 
Henderson and Vedral.) In the paper-2., the min- 
imal difference I — I' = Q has been identified 
with an amount of quantum correlation and has 
been called the quantum discord — a measure 
of the quantum excess of correlations, a mea- 
sure of the quantumness of correlations. Taking 
into account that 



I = C 



(12) 



we see the equivalentness of results of 
Henderson- Vedral and Ollivier-Zurek. 

Quantum discord displays a number of prop- 
erties (see, e. g., the reviewsi^-®). We note be- 
tween them the following ones. For a pure state, 
discord coincides with the entanglement E. In 
mixed states, the quantum correlation (discord) 
can present even in that case when the entan- 
glement is absent. Quantum discord Q > 0. 
Discord is limited from above by the entropy of 
one subsystem, Q < S'(p^(B))- In particular, if 
the system is two-qubit and the logarithm base 
for entropy equals two then Q < 1. 



III. HAMILTONIAN AND DENSITY 
MATRIX 

Consider a system consisting of two identical 
particles with the spins 1/2 which couple be- 
tween themselves by the magnetic dipolc-dipole 
interaction. Let moreover, the external homo- 
geneous magnetic field with induction B was 
applied to the system. Then the Hamiltonian 
of a model can be written as (see, e. g.^) 



n = ndd + nz, 

where the dipolar part is 



(13) 



Held — 



and the Zeeman energy equals 
1 



[(Ti(T2 - 3(n • o-i)(n • (T2)] (14) 



In these equations, fiQ is the magnetic constant 
(magnetic permeability of free space), 7 is the 
gyromagnetic ratio, (Ti,2 are the vectors of Pauli 
matrices at the sites 1 and 2, r is the distance 
between the spins in a dimer, n is the unit vec- 
tor in the direction from one spin to the other, 
and B is the vector of magnetic field induction. 
The dipole-dipole interaction reflects the exact 
law (in that sense that it does not contain the 
fitting parameters), the interaction is sharply 
anisotropic and long-acting (in contrast, say, to 
the exchange interaction). 

In the spherical coordinates when B = 
{0,0, B) and n = (sin0, 0, cos0), the Hamilto- 
nian - (|15l) takes the form 



H : 



47r Ar^ 

1 



sin t 



- 0-: 



X ((72 sin 6 - 



- o"! sin ( 



' sin 9) 



(16) 



It has been shown~ that when the polar angle 
= ^ (the external field is applied perpendicu- 
larly to the direction of dimer longitudinal axis 
the entanglement in the system is maximal. On 
the contrary, when = or tt (the longitudinal 
dimer axis is parallel to the external field) the 
entanglement between spins is absent (below we 
prove this strongly /exactly). 

Because the special interest is to discover and 
investigate the quantum correlations without 
entanglement, we will consider from this point 
the case = 0. By this, the Hamiltonian ([Tl 
takes the form 



where the dipolar coupling constant equals 
47r 2r^ 

and the normalized external field is 
h = jhB. 



(17) 
(18) 

(19) 



-Hz = -^lh{ai+a2)B. (15) 



In dipole-dipole coupled dimer, the anisotropy 
parameter is A = —2. However, below we will, 
in some cases for the sake of generality, extend 
the values of A. But all graphical and nu- 
merical material in our paper is presented for 
A = -2. 

The Hamiltonian (fT7|) corresponds to the 
XXZ model in Z field. In the matrix forn, it 
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IS given as 

A 



H 



\ 



2 ^ 

D 



D 



-D 



W + hJ 



(20) 

The 2x2 subblock presented here is a cen- 
trosymmetric matrix which under the orthog- 
onal transformation 



O 



1 

71 



1 1 

1 -1 



(21) 



undegoes in diagonal form. 

The energy spectrum of the Hamiltonian ((20] 
with condition A = — 2 consists of two indepen 
dent on external field levels 



2D, 



(22) 
(23) 



El 

and two levels 

£'3,4 = -D±h. 

Because D > 0, the ground state energy is 

Eo = -D+\h\. (24) 

In absence of external field, the ground state is 
two-fold degenerate. 

We will consider the dimer in a termal equi- 
librium state. In this case, its density matrix 
p = pab has the Gibbs form 



1 



p = — cxp(-/3-H), 



(25) 



where /? = l/kgT, ks is Boltzmann's constant, 
and Z is the partition function. 



Z = Trexp(-/3-H). 



(26) 



Performing necessary calculations we find that 
in the original (standard) basis |00), |01), |10) 
1 11) the density matrix has the structure 



/a 



P = 



V 



b V 
V b 



(27) 



where 



a= ^exp[-/3(|i?~/i)], 
d=|exph/3(fZ? + /i)], 



(28) 



and partition function is 



Z = 2(ch/?D 



e 



-0DA 



ch I3h} 



JDA/2 



(29) 



Expressions ((28)) satisfy the condition 

a + d + 2b^l, (30) 

which provides the normalization Trp = 1. 

Expanding the density matrix (|27p into pow- 
ers of Pauli matrices we obtain it in the (nor- 
mal) Bloch form 

P=\[l + ia- d){at + al) + 2v{afcT- + diVf) 



+ (1 - 4b)<jla^, 



(31) 



Expansion coefficients are the unary and binary 
correlators, 



(32) 

= l-4b=l-^e^''^/\h(3D, 

Zj 



•\"2. 



(33) 

Z ^ 
(34) 

Brackets denote the trace operation for the 
expression in brackets with density operator, 
(•) = Tr(-yo). The coefficient a — d = m in the 
expansion ((3ip is equal to the z components 
(projections) of the Bloch vectors for the re- 
duced density matrices of subsystems A and B. 
Moreover, the quantities m, G\\ , and Gj_ have a 
physical sense, correspondingly, as the normal- 
ized magnetization, longitudinal and transverse 
components of correlation matrix. Relations 



2m + G||), fe=i(l-G||), 



4 ^ 



Gil 



(35) 



give the expressions for the matrix elements of 
the density operator ((27|) through the system 
correlators. 

From Eq. ((321) -((Ml) ^e find the high temper- 
ature behavior for the magnetization and cor- 
relation functions, 



1 h 



A h f D Y 
2 ksT^^D \ Wr) 



G\\{T,h) = 

+o(i/r3). 



A D 



2 knT 



0(1/T3), 
(36) 

2' 



D 



knT 



(37) 



5 



G^iT,h) = -- 



1 D 



1 rl 
-3 



2 ksT 

2 



D 



D 



O 



1 

y4 



(38) 



Thus, when T — > cxd, the main contributions 
to the correlators tend to zero according to the 
law \/T. By this, a weak external field does 
not exercise an influence on the G|| and Gj^. 
Moreover, the leading term in the expansion of 
the transverse correlator G± does not depend 
on the anisotropy A. 

At lower temperature, the correlation func- 
tions for the dipolar dimer (r'>0,A = — 2)in 
absence of external field behave as 



G, 



1 — exp( 



keT' 



1 , D 



(39) 



(40) 



If the field h 0, the additional statistical 
weights at the exponents arise (due to a change 
of ground state of the system; see below). 



G|||t^o « 1 - 2exp(- 



Gi 



D 



D 



), (41) 



(42) 



Thus, the lower temperature behavior of cor- 
relators is described by a function of the form 



IV. QUANTUM ENTANGLEMENT 

The entanglement through the relation 



1 + - G% /I + - G2 

E = 7, log2 



1 - Vl -G2 , /l- Vl-G2^ 
log2 



(43) 



is expressed via the concurrence G^^'^^ . 

In the case of density matrix having the 
block-diagonal form 



Xi w 
W* X2 



(44) 



the equatio] 
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G = 2 max{|w| - yiZItl^, 0} (45) 

serves for calculation of concurrence. The 
expression (|45|) is a particular case of Hill- 
Wootters formula^i^ which allows to calculate 
the pairwise concurrence in general two-qubit 
system. 

Our density matrix p7| has the form of 
Eq. (jHl). Using Eqs. ^ we find that 



\w\ — y/uiU2 = —V — v ad < 0, (46) 

if A < —1. This inequality is valid for arbi- 
trary external field B. Hence, the concurrence 
(|45)) and together with it the entanglement (|43|) 
are equal identically to zero. So, in the dipole- 
dipole (A — —2) dimer under question, the 
quantum entanglement is absent for all tem- 
peratures and arbitrary longitudinal fields. 



V. INFORMATION CORRELATIONS 

In this section, we give a calculation of in- 
formation correlations in dipolar dimer both in 
and out magnetic field. 

A. Arbitrary external field 

Discord in a system with the density matrix 
having the structure ((27|) equals^i^ 



where 



Q = min{Qi,(32}, 



Qi ^ Sa~ Sab - a log2 



(47) 



a + b 



-6l0g2 ( — ^ ) - fol0g2 ^ 



-dlog2 



a + b 
d 

b + d 



b + d 



(48) 



U2 



Q2 = Sa- Sab - Si log2(5i - S2 log2(52, (49) 
and 

(5i,2 = i[l±((a-rf)2+4«2)i/2]. (50) 

In equations (gH]) and (gS]), 

Sa = -(a + b) log2(a + b) - {b + d) \og^{b + d) 

(51) 
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is the von Neumann entropy of reduced density 
matrix pA and 

Sab ~ — alog2a ~ d\og2d — {b + v) log2(fe + v) 
-(5-t;)log2(6-«) (52) 

equals the von Neumann entropy for the density 
matrix p of full system. 

The total correlation is / = 2Sa — Sab, and 
classical one C = I — Q. These relations and 
also (gZll-dSll) together with Eqs. ^ and ^ 
define the quantum discord, classical, and to- 
tal correlations as functions of the temperature 
and external magnetic field by arbitrary values 
of parameters D and A. 

At high temperatures. 



Qi 



1 



D 



41n2 \ kBT 



A 



D 



81n2 \ kBT 



-o(i/r4), 

(53) 



Q2 



81n2 



(1+A2) 



D 



-Oil/T^). (54) 



Thus, the main terms of these expressions do 
not depend on the external magnetic field. 
From the relations ([5^ and ([5^ . one can see 
that by |A| > 1 the discord, according to 
Eq. (|T7)) . is defined by the branch Qi. Hence, 
when A>lorA<— 1, the quantum discord 
by high temperatures behaves as 



QI 



1 



1 



41n2 [kBT/Df 



(55) 



Thus, it does not depend on both external mag- 
netic field or interaction anisotropy A. With 
increasing the temperature, the quantum cor- 
relations decrease according to the law 
i. e., essentially rapidly than ordinary statisti- 
cal correlations which, how it can was empha- 
sized above, tend to zero according to the law 
l/T. 

From equations presented above, we estab- 
lish also that by high temperatures the classical 
correlation is 



1 



(56) 



%\n2 [kBT/Df 
and does not depend on external field. 

B. Dimer in absence of external field 

In important particular case of zero external 
field, the information correlations can be calcu- 
lated via the simpler formulas'^-. In absence of 




kjT/D 



Figure 1: Quantum discord in dipolar system 
(A = —2) as a function of temperature and ex- 
ternal longitudinal magnetic field. 



field, a — d and Sa — Sb = 1- Besides, using 
expressions (155)1 and it is not difficult to 

establish that by A < — 1 one has Gn > |Gj_| 
(below this will be seen in graphics) . Therefore, 
in absence of magnetic field, the classical part 
of mutual correlations is 

C = l [(l+G|||)log2(l+G|||)+(l-G,i|)log2(l-G,i|), 



and the quantum discord equals 
2G 



(57) 



g = - [(I + 2Gi - G||)log2(l + 2Gi - G||) 

-2(1 - G||)log2(l - G||) + (1 - 2G^ - G||) • 
log2(l-2Gi-G||)]. (58) 



How can one see from Eq. ([57)) . the classical 
correlation is completely determined by the lon- 
gitudinal correlator Gy. Using expressions for 
the correlation functions ([55)) . (15^ and setting 
/i = in them, we get the following formula for 
the discord, 



(ch(^)+exp(-i^))ln(2) 



(59) 



For a dipolar dimer [D > 0, A = —2) in zero 
field and at lower temperatures, we have 



1 



exp{-D/kBT). 



(60) 



Arbitrary order derivatives with respect to the 
temperature for the function in the right hand 
of this equation is zero at T = 0. Therefore, by 
small deviation of temperature from absolute 
zero, Q « 0. This is connected with existence 
of a gap in the energy spectrum of the system. 
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5 1 1.5 2 2.5 

kT 
D 



Figure 2; Isolines of quantum discord in dipolar 
system (A = —2). 



VI. DISCUSSION 

A behaviour of quantum discord in the 
dipole-dipole system under question (A = —2) 
is shown in Fig. [TJ We see the smooth hill-hke 
surface stretched in the temperature axis direc- 
tion. Along a straight hne h = 0, two ridges 
go; one ridge goes in the direction of tempera- 
ture decreasing and the other in the direction 
of its increasing. The surface is symmetrical 
under reflection in vertical plane going through 
the hill top and the straight line ft. = in the 
temperature-field plane. By given temperature, 
the discord is maximal in absence of external 
magnetic field, that is, the field leads only to a 
suppression of quantum correlation. At the ab- 
solute zero temperature, the discord in dimer is 
identically equal to zero. In the high tempera- 
ture limit, quantum correlation is also vanished. 
In Fig. 2, the cross sections (profiles) of a dis- 
cord surface are shown by different values of Q. 
It is seen that the isolines form a set of non- 
crossing ovals. 

To study in detail the behavior of discord, we 
represented in IVII its temperature dependence 
by different values of external field. The discord 
reaches the largest value at /i = and 

fcsT^/D^ 0.881297... . (61) 

The discord in this point equals 

Qm = 0.083 061... . (62) 

what is 8.3% of maximal value which is possi- 
ble in any two-qubit system. The value ()6ip , we 
have found both by a numerical search of max- 
imum for the function Q of kBT/D and from 



Q 









A = -2 


- 






- 






/ \ 

2 
















/ 













0.5 1 1.5 2 2.5 3 

kaT/D 

Figure 3: Discord isoterms by different values of 
external field. 



solution of transcendental equation 

^(e'^'^-l-chx-l-A-sha;) = (shx + A-cha;) In(chx), 

(63) 

where x = D/kBT. This equation follows from 
the condition dQ / dT = in which the function 
Q{T) is given by expression (f59l) . 

One can see from Fig. 3 that, with increasing 
the field, the discord maximum is moved in the 
direction of higher temperatures. This is a good 
feature. However by this the value of discord at 
the maximum is less than its value at the same 
temperature in absence of external field. 

We see also that with the help of external 
field and temperature (these parameters are in 
our hands), one can control the discord varying 
its value from zero to 0.083 bit per a dimer. 

Let us consider the case when the dipolar 
dimer is in absence of external field. The tem- 
perature dependences of ordinary spin-spin cor- 
relation functions Gy and {G± < 0, there- 
fore its absolute value is taken), as well as total 
classical correlation C and discord Q are pre- 
sented in Fig. 4. At zero temperature, both 
longitudinal and total classical correlation are 
equal to their maximal magnitudes. The tem- 
perature increasing leads only to their decrease. 
Both the value of transverse correlation |G_l| 
and quantum discord are zero at T = 0. 

To clarify the situation in correlation behav- 
ior at zero, we consider the limit of density ma- 
trix (P7|) when r — 7> 0. Using the expressions 
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Figure 4: Temperature behavior of statistical (Gy 
and |Gx|), classical (C), and quantum (Q) correla- 
tions in absence of magnetic field. 
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Figure 5: The same as in Fig. 4 but in presence of 
a field. 



for its matrix elenients ([28|) one finds 

p\t=o \ j . (64) 

This matrix has the expansion 

P= ^(|0)(0U®|0)(0|b + |1)(1U®|1)(1|b, (65) 

i. e., it is a sum of the direct products of sepa- 
rate particles. Such state belongs to the class of 
purely classical ones^ and therefore, according 
to the criteria'^^, all quantum correlations in it 
are entirely absent. This explains that we ob- 
tained Q\t=o — 0. On the other hand, classical 
correlations are available and equal the total 
correlations C = I. 

So the behavior of system at T = is purely 
classical. However, with increasing the tem- 
perature, the behavior acquires the quantum 
features, i. e., the temperature leads to the 
generation of quantum correlations. Such un- 
usual phenomenon (it is much ordinary when 
the temperature destroys the quantum states) 
one can explain by following. In the Hamil- 
tonian p7)) with A < — 1, the purely classical 
ferromagnetic contribution 

H,^ = ADafa^ (66) 

is dominate at lower temperatures. The trans- 
verse components of spins are actually "frozen" 



(the value G± — witnesses this). With in- 
creasing the temperature these degrees of free- 
dom "revive" and the system from classical be- 
comes a quantum one. But when the temper- 
ature fluctuations begin to exceed the system 
energy gap, the temperature produces its usual 
destroying action and the correlations go down 
to non. 

From Fig. 4 and asymptotical behavior of 
correlations, it is not difficult to notice a corre- 
spondence in qualitative behavior, on the one 
hand, C and and, on the other hand, Q and 

G5_- The latter correspondence is not an ac- 
cident. Indeed, if the quantum correlation is 
measured by geometrical discord Qg, i. e., by a 
distance (in sense of the Hilbert-Schmidt norm) 
from the given state p to the nearest classical 
one, the general formula^i^ applied to the den- 
sity matrix pip by an absence of external field 
yields 

Qg = Gl. (67) 

Consider now the situation in a presence of 
external magnetic field. In Fig. 5, the behavior 
of different correlations is again shown but for 
h/D = {).!. It is seen that correlations G||, G±, 
and Q were changed weakly. But the behavior 
of classical correlation G undergoes an essential 
change — it tends now to zero when T — ^ 0. At 
the point of absolute zero temperature, there 
are no any correlations in the system at all. 

To understand a happening, we turn out 
again to the density matrix. The matrix (|27p 
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Thus, any statistical relation between spins is 
also absent. 

Fig. 6 shows the behavior of classical corre- 
lation C upon the external field. The curves 
have the form of bell-like splashes. Their largest 
value is arrived at in the point when the field 
vanishes. With decreasing the temperature, the 
maximums become more narrow and their value 
tends to the value C = 1. At T = 0, a splash 
becomes infinitely thin. 



(72) 



h/D 




Figure 6: Classical correlation vesus the field by 
difi'rrent values of temperature. 



in the limit T — )■ is equal now to both 
/I 



p\t=o 











(68) 
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if /i > or 



p\t=o = 



V 



(69) 



1 



if ft, < 0. These states correspond to the com- 
pletely orderings of spins along the field. Both 
matrices and (p^ are factorized into a di- 
rect product of two density matrices 2x2, 

/I 

















I) 



V 



(70) 



1 



This means that the states (|68p and ([S^ are 
completely uncorrelated and therefore any cor- 
relations are absent in them. 

Emphasize that Gy 7^ at T = does not 
contradict to the said above. In the factorized 
state must be zero only the centered correlator. 
By h 0, the spins at the point of absolute 
zero temperature are ordered and therefore 

(K - K»(^2 - (^l»> = 0. (71) 



The external field sweeping will on a moment 
lead to a spasmodic jump appearance of classi- 
cal correlation and at the same time to a disap- 
pearance of this correlation by going the field 
of the zero point. By this, the quantum corre- 
lation in the system does not arise. 

Let's turn to the available experimental data. 
Measurements performed by the nuclear mag- 
netic resonance (NMR) at room temperature 
show that in gypsum crystals CaS04 • 2H2O 
the distance between protons in each water 
molecule is r = 0.158 nm^-'^ (see also, e. g., 
the book^'^). For protons, the gyromagnetic 
ratio, as it is known^^, equals 7 = 2.675 • 
10® rad/(c-T), therefore the dipole-dipole cou- 
pling constant (jl8p in gypsum is D/ks = 
0.73 /iK (in temperature units). Consequently, 
in accord with (j6ip , the maximum discord value 
Qmax = 0.083 will arrive at the temperature 
0.64 /iK. At room temperature (T = 300), 
the discord in gypsum, according to ([55]) . must 
be equal to Q ~ 2 • 10~^®. In spite of ex- 
tremely small value of quantum correlations in 
spin-nuclear systems at room temperatures, at 
present the attempts are undertaken to detect 
their by NMR methods^^ii. 

As an another example, let us consider 1,2- 
dichloroethane CIH2C — CII2CI. In this com- 
pound, two protons at each carbon atoms are 
coupled much stronger between themselves by 
dipole-dipole interaction than with protons be- 
long to an other carbon atom. NMR mea- 
surements performed on solid dichloroethane at 
the temperature 90 K have shown that here 
r = 0.17(2) nm.'^-^ (see also^i^). Using again 
the relations (|18p and (pT|) . we estimate the 
temperature for the discord maximum in this 
substance as T„i — 0.517 //K. At the temper- 
ature 90 K, the value of quantum correlations 
must equal Q ~ 1-5 • 10" 
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VII. CONCLUSIONS 

In the paper, a study of information correla- 
tions in dipolar dimers both in absence of exter- 
nal magnetic field and in the field directed along 
the longitudinal axis of a dimer has been per- 
formed. It bas been shown that the quantum 
correlations are completely absent at the ab- 
solute zero temperature and arbitrary strength 
of magnetic field but they arise with increasing 
the temperature. We proposed an interpreta- 
tion such of phenomenon. 

In the high temperature region, the discord 
obeys the law 

Q^Gl^ l/{kBT/Df. 

This, in particular, means that at room, ni- 
trogen, or helium temperatures, the non-zero 
quantities of spin-spin correlations between the 
XX or yy spin components can serve as a witness 
of quantum correlations. The spin-spin correla- 
tions go down with increasing the temperature 
according to the essentially slow law and 
one can directly measure them, for example, in 
scattering experiments. 



It was shown also that the classical correla- 
tions have a sharp maximum in the point of 
zero external magnetic field when T — > 0. 

A qualitative analogy between the quantum 
discord and the squared correlator Gj, has been 
found. In supporting this observation, it was 
shown that the geometrical discord equals Qg = 
Gi. 

In the paper, the estimates for the thermal 
quantum discord between spins of hydrogen nu- 
cleous in gypsum and 1,2-dichloroethane 
have been done. 

Our investigation can be extended to many- 
nuclear clusters if to perform a density matrix 
reduction for all spins except any two. Exper- 
imental NMR data are available for clusters in 
a form, for example, triangles, tetrahedron, lin- 
ear magnetic structures, etc.j^i^i^. 
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